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We develop the Ginzburg-Landau theory of the vortex lattice in clean isotropic three-dimensional 
superconductors at large Maki parameter, when inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov 
state is favored. We show that diamagnetic superfluid currents mainly come from paramagnetic 
interaction of electron spins with local magnetic field, and not from kinetic energy response to the 
external field as usual. We find that the stable vortex lattice keeps its triangular structure as in 
usual Abrikosov mixed state, while the internal magnetic field acquires components perpendicular 
to applied magnetic field. Experimental possibilities related to this prediction are discussed. 

PACS numbers: 74.20.Fg, 74.25. Dw, 74.62.Dh 



I. INTRODUCTION 

Orbital and paramagnetic effects are both important in 
the suppression of superconducting state. While orbital 
effect leads to the formation of Abrikosov vortex lattice 
below the orbital upper critical field H c2 o « <I>o/27r£o 
in type II superconductors,— paramagnetic effect deter- 
mines the paramagnetic limit of superconductivity H p = 

A /V2^r^ an d promotes the tendency to the Cooper 
pairing with nonzero momentum - so called Fulde-Ferrell- 
Larkin-Ovchinnikov state (FFLO)4£ Here $o is the flux 
quantum, £o is the superconducting coherence length, 
A is the superconducting gap at zero temperature and 
fi = g(iB/2 is the electron magnetic moment. 

The interplay of both effects takes place at large 
enough ratios \/2H C 2o/H p = a Mr called the Maki 
parameter^ In this regime, the (H, T)-phase diagram 
of isotropic s-wave superconductor in the clean limit 
was studied by means of a Ginzburg-Landau (GL) 
functional^ Critical assumption was made that screen- 
ing supercurrcnts are not important, and the local field 
in the superconductor was taken equal to the external 
field. This is true when the GL parameter k is large 
enough. For the clean superconductor, the GL parame- 
ter is related to the Maki parameter 



^k F r e (m* /m) ' 
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here fcj? is the Fermi momentum, r e = e 2 /mc 2 is classical 
radius of electron and m* /m is the ratio of an effective to 
the bare electron mass. It is clear from this relation that, 
if the Maki parameter is large, the GL parameter has even 
larger value, hence the assumption done in Ref. H seems 
reasonable. 

If, however, one is interested in the magnetic response 
of superconductors in FFLO state, it is necessary to re- 
lax this assumption. In the present article, we develop 
the theory of FFLO state in a clean isotropic three- 
dimensional superconductor, including the space varia- 
tions of currents and fields. 

The article is organized as follows. In Section II, we in- 
troduce the free energy functional including the Zceman 



interaction with spatially non uniform magnetic field. 
Being unimportant for the upper critical field determi- 
nation (Section III) and for the vortex lattice structure 
found in the Section IV, this term is crucial for the cur- 
rent and field distribution in FFLO modulated supercon- 
ducting state (Section V). There we find that the internal 
magnetic field has components perpendicular to applied 
magnetic field. The experimental possibilities related to 
this theoretical prediction are discussed in Conclusion. 



II. FREE ENERGY 



We investigate the phase diagram of a superconductor 
near the upper critical field, but far from critical temper- 
ature T c . For this purpose, a GL expansion of the free 
energy in powers of the order parameter and its gradients 
is possible provided the length scale of variations of the 
order parameter, determined by the magnetic length, re- 
mains large compared to the superconducting coherence 
length. This is indeed the case in the limit of large Maki 
parameter when the critical field is mainly determined 
by paramagnetic depairing effect. 

As in Ref. ©, the free energy density was derived 
both in frame of Gor'koviS and Eilcnbcrgcr— or Larkin- 
Ovchinniko v 12 ! 13 formalisms^ 
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Here, T n § is the free energy density in normal state in 
absence of magnetic field, D = — iV + 2eA, h = rotA is 
the local internal magnetic field and the coefficients in the 
functional depend both on temperature T and induction 
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B determined by the spatial average h = B = Bz: 

fl-^., «~^*., 

_*£ K e = -»/V„pLj, (3) 
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where Vo = ?n*fci?/27r 2 is the normal density of states at 
Fermi level, vp is a Fermi velocity, to = ixT(2v + 1) is a 
Matsubara frequency, and 

K n = 2T9? V — , L n = 2T3 V - — 

(4) 

Standard form of GL functional is given by the terms 
in the first line of eq. ([2]) only. In the paramagnetic limit, 
transition from normal to uniform superconducting state 
takes place at critical field B C (T) defined by a = 0. Along 
this transition line, the coefficients (3 and 7 become nega- 
tive at T < T* ~ 0.56T C . This signals possible instability 
toward FFLO state with spatial modulation of the order 
parameter A, as well as possible change of the normal to 
superconducting phase transition order, at the tricritical 
point (T*,B* = B C (T*)). Higher order terms must be 
retained in the functional density © to consider these 
effects, as shown in Ref. 0j| 

A functional similar to ([2]), but including arbitrary dis- 
order, shape of the Fermi surface and pairing symmetry 
of the superconducting state was derived in the limit of 
k — > ex).- At finite value of K, the coordinate dependent 
deviation of magnetic field from external field manifests 
itself not only in gradient terms, but also in Zeeman in- 
teraction with electron spins. This results in the last 
term in the functional which is absent in Rcfs. 

mas 

and which simply corresponds to local decrease (enhance- 
ment) of the critical temperature T C (B) when h z (x) > B 
(< B). Note that the coefficient e is proportional to B. 
Hence, the corresponding term in functional @ is negli- 
gibly small in ordinary GL region near T C (B — ► 0) 

III. UPPER CRITICAL FIELD 

At the second order phase transition between normal 
and superconducting state, the magnetic field is uniform, 
h C 2 = Bz. The linearized gap equation obtained from (|2|) 

qA + 7 D 2 A + 5[(D 2 ) 2 + A- 2 ]A = 0, (5) 

where A -1 = y2eB is inverse magnetic length, is solved 
with: 

A = i p (x,y)f(z), D 2 A=(± + q^jA (6) 

where <fo(x, y) is the linear combination of Landau wave 
functions with level n = for the particle with charge 2e 



under magnetic field B, multiplied by exponentially mod- 
ulated function along z-direction, f(z) = e ±lqz . (Note 
that, in principle, higher Landau levels could also be con- 
sidered. But they may be realized only at low temper- 
atures when transition from normal to superconducting 
state with lowest Landau level has turned first order—, 
see Section IV.) 

At large Maki parameter, the upper critical field is 
close to B C (T). We note that a ~ e{B - B C (T)) and 
eq. ([5]) then defines a field: 

B(q) = B c (T)-l(\- 2 + q 2 )--{(\- 2 + q 2 ) 2 + \-% (7) 

e e 

The upper critical field B c i is found by taking the maxi- 
mum of B(q) with respect to q. Analysis of eq. {7J shows 
that, in the presence of orbital effect, the critical temper- 
ature T* below which FFLO modulation appears, defined 
by 7 + 25/ A 2 = 0, is decreased compared to its value T* 
in absence of orbital effect £ 

f*~T* -1.2—. (8) 

OiM 

Namely, at T > T*, the usual superconducting state with 
q = appears with critical field 

Bc2iT) = Bc(T) -tEfh^j^m. (9) 

While, at T < T*, the FFLO state appears with finite q: 
-f + 25(q 2 + \-' 2 ) = 0, (10) 

and critical field 

IV. PHASE DIAGRAM AND VORTEX 
LATTICE STRUCTURE 

In order to determine the structure of the vortex 
lattice state, following variational procedure similar to 
Refs. HH, consideration of higher order terms in free 
energy density {2j is required. Just below the upper 
critical line defined by B c2 (T), the magnetic field is 
partially screened by supercurrents and we decompose 
h = B + hi, with B < B C 2 and hi = 0, and, corre- 
spondingly A = A + A\. At T > T*, conventional 
Abrikosov state (A state) is realized, and f(z) = l. At 
T < T* , FFLO state is realized and two possible modula- 
tions could appear at B < B c2 (T): the so-called FF state 
with exponential modulation, f(z) = exp(zgz), and the 
LO state with sinusoidal modulation, f(z) = v^singz. 

By minimizing free energy, we may determine the lat- 
tice geometry^ and, below T*, the temperature range 
when FF or LO state is favored. For this, we write spa- 
tial average of the free energy density ( [5]) in the form: 

Tl = B ' ' + Hl + MA, A) + J%E, A) (12) 
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where Ti and T± collect together quadratic and quartic 
terms with respect to A, respectively. By making substi- 
tution A — * (1 + e)A in T s and requiring that linear-in-e 
terms vanish^ we get: 



= T 2 {A,A)+2T A {A 1 A) 



(13) 



^ 2 (A,A ) 



5T 2 



A 1 .^(A ! A ) + 2^ 4 (A,A ). 



By variation of the free energy with respect to A, we get 
the Maxwell equation relating internal magnetic field to 
screening currents: 

i-rotA 1= j s = -^(A,A ), (14) 

up to second order terms A. We insert Eq. (Tl4|) into (Tl3|) 
and we integrate the second term by part. Then, eq. (fH 
yields: 



IAP 



jr 2 (A,A )/|A|2 



^4(A,A )-h?/87r /(|A| 2 ) 



(15) 



where r.h.s depends on the structure of order parameter 
only. Near transition, we note that JF 2 (A, Ao) ~ e(B — 
B c2 (T))|A| 2 . Inserting eq. (Jl5j) into ([12]), we obtain: 

e 2 (B - B c2 (T)) 2 



— B 2 



(16) 



4[^ 4 (A,A )-^/87r]/(|AP)2 

The equilibrium vortex lattice structure is thus the one 
which minimizes the denominator of the second term in 
r.h.s. of eq. (fl6|) . 

At large enough GL parameter, the denominator in the 
second term of eq. (|16[) is dominated by its first term. 
Noting that the gap function ([6]) obeys the properties: 



(-D±A) 2 
and using eq. 

■F 4 (A,A ) 

(W) 2 

where 
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TUl) in FF or LO state, we find: 
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in A state, 

in FF state, 

in LO state. 

(18) 

(19) 



is the Abrikosov parameter.— Comparison of the expres- 
sions in r.h.s. of (TT5|) yields^ that FF state is realized in 
the temperature range defined by 9/28 > \ 2 K 3 jv 2 F K^ > 
3/10, that is, T = f* - 0.08T c /a M < T < f*. LO 
state is realized below these temperatures, and the tran- 
sition from normal to LO state changes its order at 
X 2 K 3 /v 2 F K 5 > 3/4, that is T < T x = f* - 2T c /a M . 

At all temperatures, free energy (|16p is minimized 
when (3a is minimal, that is for triangular vortex lat- 
tice (with (3a = 1.1596). Corresponding phase diagram 
is shown qualitatively in the Figure. 
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FIG. 1: Qualitative phase diagram of a clean three- 
dimensional superconductor with large Maki parameter. Thin 
(thick) line is for second (first) order transition. Transitions 
shown with dashed lines are not discussed in the present work. 



V. ELECTRODYNAMICS OF VORTEX STATES 

To consider the situation at finite GL parameter we 
need to evaluate the term h\ in eq. (| 16[) where hi solves 
the Maxwell equation (fT4|) with: 



Jkir 
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-|rotrot[A(DA)*] + c.c. } + 8e 2 5B x V|A| 2 , 



Jz 



rot(|A| 2 z). 



(20c) 



Here, j'km originates from usual superconducting kinetic 
energy response to the spatially varying magnetic field. 
The Zeeman current jz, arises from the interaction en- 
ergy of the superconducting diamagnctic correction to 
the normal- metal paramagnetic moment, with the spa- 
tially varying magnetic field hi. 

Making use of the form ^ for the gap and its property 
(at Landau level n = 0) 



A*(JD ± A) + c.c. =rot(|A| 2 z), 



we find: 
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-rot 
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where e = e + 2e(7 + 4<j/A 2 ) in A state, and, with use of 
eq. (TTU]), e = e + Ae8/X 2 in FF and LO states. And for 
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longitudinal component, again with use of eq. (|10[) , we 
have in the lowest order in |A| 2 : 



3s,\\ 



in A and LO states, , . 

-4e^V 2 |A| 2 z in FF state. { ' 



Let us note that in A state or LO state, the supercur- 
rents flow in planes perpendicular to the external field 
only, while in FF state, they also flow in parallel direc- 
tion. Close to the critical temperature T*, the simple 
estimation shows that 



Jkij 



Jz 



1 



a 



(23) 
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Hence, we may neglect jkm in perpendicular component 
of current. In A and LO states, the Maxwell equation 
(fT4|) acquires the following form 



rot/ii 

47T 



e rotlAI 2 ; 



(24) 



In FF state, we may keep the small term corresponding 
to the parallel-current component and we get: 



rot/ii 

47T 



rotlAI 



^V 2 |A| 2 z. (25) 



Let us remind that in conventional Abrikosov state in 
vicinity of critical temperature, where coefficient 7 is pos- 
itive, the supercurrents are diamagnetic: that is, they 
create a magnetic moment directed opposite to the ex- 
ternal field. The value of e is always positive, hence as 
in the Abrikosov case the orbital currents in FF and LO 
modulated states are also diamagnetic despite their Zee- 
man origin. 

Let us now determine the distribution of fields in vor- 
tex state. The component h\ of the magnetic field is 
periodic and is found from Maxwell equations (|24[) or 
(pjj) . and 



div/ii = 



(26) 



together with the condition h\ = 0. In A state, these 
equations are solved with: 



-^=e(|A| 2 -|A| 2 )z. 
In FF state, they are solved with 



(27) 



hi 



4e5q 



- -L = £ (|A| 2 - |A| 2 )i + -y* (z x V|A| 2 ). (28) 
In LO state, we search for the solution in following form: 



- ^ = e(|A| 2 - |A| 2 )z + 2qV[ X (x, y) sin2gz], (29) 

and we find that x( x tU) is an auxiliary function which 
solves: 



Vlx + 4g 2 x = e|(/3 | 



(30) 



In FF and LO states, in contrast with the situation in 
Abrikosov vortex state, the internal field has components 
perpendicular to applied magnetic field. 

Evaluating h\ in A and FF state, we get (neglecting 
the transverse component of field in FF state) : 



h 2 



8tt(|A| 2 ) 2 
In LO state, we get: 
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(31) 



(32) 



A rough evaluation of the last term in (f3"2"|) can be done 
by estimation of — V\x ~ ^ 2 X i n (|30|) . Consequently 
similar to (f3"Tj) we obtain^: 
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2ne 2 (Cp A - 1), 
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where in A and FF states C = 1 and for LO state 
3 2q 2 \ 2 



Clo 



2 l + 4g 2 A 2 ' 



(33) 



(34) 



Inserting the result for hi in (fTB]) and using eq. (Ti"%j) . 
we can present the free energy in usual form>iii£ 



(B-B c2 (T)) 2 



8ir 87r[l + (3 A {2 K 2 cff -C)Y 



(35) 



with effective, temperature-dependent, GL parameter 
K ff defined by 



iVo 
4e 2 



fK 4- in A state, 
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(36) 



We note that the form (|35|) of the free energy requires 
> 0, that is T > T x - Below T x , normal to FFLO 



l cff 



state transition becomes of the first order and requires 
higher order terms in the gap to be retained in eq. ([2]). 
Moreover, at T > T\, in analogy with type I/type II su- 
perconductors, the free energy (|35|) is indeed minimized 



with vortex lattice state only if «; e ff > y/C/2. The situa- 
tion in vicinity of the point «„ff = y/C/2 requires special 
investigation similar to Ref. [l8| at T — > T c . Evaluating 
Ateff at temperature T* corresponding to the tricritical 
point for N/A/FF states, we find: 



\rff ■ 



(37) 



A[ 



The theory of vortex lattice in FFLO state thus applies 
at: 
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> 1. 



(38) 
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We now determine the diamagnctic response. Applied 
magnetic field is found from thermodynamic relation: 



H = Ait 



dB ' 



(39) 



From relation B = H + AttM, we obtain the magnetiza- 
tion induced in superconducting state at a given applied 
magnetic field 



M(H) = - 



B, 



c2 



H 



4tt (2^ ff _ C)(3 A 



(40) 



We note that the derivative of induced magnetization in 
superconducting state with respect to applied magnetic 
field varies with the temperature along the upper critical 
line (T, H C 2(T)). Its most peculiar features are: I) it 
increases abruptly close to temperature To, as C = 1 at 
T > Tq when the transition is from N to FF state, and 
C > 3/2 at T < Tq when the transition is from N to LO 
state; 2) it diverges at T < T\ when transition into LO 
state becomes of the first order. Feature 1) takes place 
close to the triple point for coexistence of N, FF and LO 
state, and it is related to the nature of FF/LO transition 
which is of the first order—. 

The determination of magnetic field distribution is also 
important for nuclear magnetic resonance (NMR). Due 
to different field distributions (f27)). ([28]) . and ([29]) in A, 
FF, and LO states, respectively, we may expect distinct 
NMR line shapes in each of theses states. 



VI. CONCLUSION 

The experimental search for FFLO state is not so easy 
due to the absence of a particular feature distinguishing 



it from the ordinary Abrikosov mixed state. We already 
pointed out peculiarities of the magnetization and NMR 
in FFLO state. More importantly, the solutions (|28p and 
([2"9"| for the field distribution demonstrate that the inter- 
nal field has components perpendicular to applied mag- 
netic field both in FF and LO states, in contrast with the 
situation in Abrikosov vortex state. The transverse com- 
ponent of the oscillating field in FF state is negligibly 
small in comparison with longitudinal oscillating com- 
ponent, while both components of oscillating field have 
comparable value in LO state. 

The effect could be revealed experimentally by means 
of small-angle scattering of neutrons polarized parallel to 
the external field. Namely, the transition to the LO (and 
possibly to FF) state should manifest itself by the strong 
increase of scattering with neutron spin flip. Another 
possibility to reveal LO and FF state is related with ap- 
plication of /iSR technique by making measurements of 
relaxation rate of precessional motion of muon spins po- 
larized along the external field direction. 

Finally, it should be noted that the effect of appear- 
ance of space oscillating transverse field component in 
LO and FF states found here for the isotropic s-wave su- 
perconductor has model independent character and will 
also be present in anisotropic materials with a different 
type of superconducting pairing. 
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